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A theory of dielectric relaxation is presented for a generalized six-site model where the 
transition probabilities for the turning of dipoles from one orientational site to another are 
arbitrary. The sites are arranged in three dimensions. The nature of the set of discrete 
relaxation times generated by this model is examined in a general manner. It is found 
that the relaxation times are all real and positive and the range of the relaxation times is 
established. 



1. Introduction 

Previous work [l] 1 with two simple models of a 
rather restricted nature has been sufficient to indi- 
cate that a set of discrete dielectric relaxation times 
can arise in a lattice consisting of dipoles that may 
occupy orientational sites arranged in three dimen- 
sions. As in studies on the single-axis rotator 
[2, 3], the source of the set of relaxation times was 
the anisotropy in the crystalline field, which hin- 
dered the turning of dipoles from one orientational 
site to another. 

In the present paper the purpose is to extend 
the treatment of the three-dimensional rotator prob- 
lem discussed in reference [1] by constructing the 
rate equations for a six-site model with arbitrary 
transition probabilities in order to examine the 
general nature of the set of relaxation times. 

2. General Six-Site Model 

2.1. Model 

A model in which a dipole of moment m on each 
lattice point can orient in any six directions is con- 
sidered. These equilibrium orientations are along 
the Cartesian coordinate axes, x, y, and z (fig. 1 (a)). 
As in previous work [1, 2, 3], the following conditions 
are assumed: 

(a) There is no correlation in the motion of dipoles 
on different lattice points. 

(b) The barrier system is the same for all molecules 
in the lattice. 

(c) The only elementary process is a reorientation 
to an adjacent site (single jump hypothesis, S-J). 
It should be noted that the general results obtained 
do not depend upon the orthogonal site orientations 
exhibited in figure 1 (a). This model was selected 
because it serves as well as any other for the general 
discussion and it was desirable to derive certain ex- 
plicit expressions useful for polarizabihty calcula- 
tions for this special case. 

i Figures in brackets indicate the literature references at the end of this paper. 



The transition probability for the turning of a 
dipole from a site i to an adjacent site j is given by 
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where ^ is a frequency factor, k is Boltzmann's 
constant, and T is the absolute temperature. The 
quantity W ti is the local free energy of activation 
barrier, which must be surmounted by a dipole in 
turning from site ^ to sitej (fig. 1 (b)). The transi- 
tion probability for the jump j to i, k jt , is given by 
eq (1) with W jt in place of W tj . As transitions are 
permitted from a site i to four adjacent sites, there- 
are 24 elementary process transition probabilities 
kij with i, j=l, 2, ... 6, and i^j in the general 
model. From physical considerations it is shown in 
the appendix (section 3.1) that of these 24 constants 
only 17 are independent, and that relations exist 
for the k tj corresponding to a closed circuit of tran- 
sitions and the inverse circuit. Thus, if we consider 
the series of site-to-site transitions 1— >2, 2— >3, 3—»l, 
and the inverse circuit 1— >3, 3^2, 2— >1, we find 2 



M 2^*23^31 — k\ 3/^32^2 



(2a) 



Other relations such as the above must apply, one 
for each site triplet for the remaining seven octants: 



^13^34^41 — '^14^43^31 
^15^52^21 — ^12^25^51 ) 



(2b) 
(2c) 
(2d) 



the four remaining equations, (2e) through (2h) are 
obtained from the equations above by substituting 
the subscript 6 for 1. Any one of the set of eight 
relations in eq (2) may be derived from the other 
seven. 3 These are the basic equations which re- 
strict the arbitrary choice of the k tj for the six-site 
model. Equations (2) are also needed for a mathe- 
matical proof noted later. 

2 A single relation of this type was given in reference 3 for single-axis rotators 
(seeeq (7) and (25)). 

3 For example, eq (2c), ^62^23^30=^63/^32^26, may be obtained as the quotient [eq 
(2a) X (2b) X (2c) X (2d) ]/[eq (2f) X (2g) X (2h)] . 
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ANGULAR POSITION IN Y-Z PLANE 

Figure 1. General six-site Cartesian axes model. 

(a) The equilibrium orientations of the dipole moment ju. (b) Local free- 
energy barrier as a function of angular position in the y-z plane through sites 

1, 2, and 6. Wn is the barrier for a transition from site 1 to site 2 and Wn is the 
barrier for the reverse transition. Vi and V2 are the free energies for sites 1 and 

2, respectively, relative to an arbitrary zero. 



2.2. Rate Equations 

Assume a lattice with a total of N dipoles where 
the number of dipoles oriented in the ith site at any 
time t is given by N t . The summation of the N t 
is always N. The number of dipoles turning per 
unit time from a given site i to an adjacent site j is 
kijNi, and the number entering site ^ fromj is k^Nj\ 
the total leaving site i to turn to adjacent sites is 
^2kijNi. Accordingly, the differential equations for 



the change in population of each site are 4 

-(h 2 +k u +h,+ku)N 1 + 

k 21 N 2 +k 3l N 3 +knN,+k 5l N 5 
= k 12 N l ~(k 2l +k 23 +k 25 +k 26 )N 2 + 
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=k 2Q N 2 +k m N 3 +h 6 N 4 + 



h 6 N 5 — (^62 + ^63 + ^64 + ^65)A 7 ' 6 



(3) 



The molecular relaxation times, rp, will be included 
in the solutions of the above equations. 

2.3. Characteristic Determinant and Its Properties 

Solutions of linear homogeneous differential equa- 
tions of the type shown in eq (3) are of the form 



ft 



(4) 



where the symbol j$ (k) represents a function of some 
or all of the k Xj . The functions e f v {k)t are referred 
to as decay functions, since the fp(k) are always 
negative for 0^2 (see below). The index ft refers 
to the mode of decay. Each decay function indi- 
cates the rate at which a particular mode of decay 
promotes the attainment of equilibrium after the 
abrupt removal of a disturbance. The fp(k) are 
obtained by forming the characteristic determinant 
of eq (3) and solving for its eigenvalues. This 
characteristic determinant is 



A(D,k u ) = 
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(5a) 



1 By the S- J hypothesis, ku=hi = k2i= . . . = /c53=0. 
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(5b) 



where we let D=d/dt and define A\, A 2 , A 3 , . 

A 1 = D+k 12 +k l3 +k u +k l5 , 
A 2 =D+k 21 +k 23 +k 25 +k 2Q , 
A 3 =D+k 31 +h 2 +k M +k m , 

AQ=D+k 62 +k 63 +k 6i +k 65 . 



The characteristic determinantal equation, obtained 
by setting A equal to zero in eq (5a), has six roots, or 
eigenvalues, fa for the operator D. The nature of 
the elements in A permits one to answer some 
questions regarding the eigenvalues and, hence, 
about the relaxation times, rp. 

The following questions may be raised: Are the 
eigenvalues real, what is their sign, what is their 
range, and are they distinct for the most general case, 
i.e., when the k i} are distinct? The eigenvalues can 
be shown to be real by applying a theorem due to 
Goldberg. 5 The line of reasoning is briefly thus: 
From the theory of matrices the characteristic 
determinant of a symmetric matrix has real roots. 
Hence, if a symmetric matrix, S, can be constructed 
that has the same characteristic determinantal 
equation as A, i.e., A=0, then the eigenvalues of A 
will be real. Sufficient conditions for forming such 
a matrix S with the same roots as A are the set of 
relations given by eq (2) and the inequalities 
kij>0. Since eq (2) and k tj >0 are imposed by 
physical considerations, it follows that the eigen- 
values are always real. 

The sign and range of the eigenvalues jp are 
specified by applying a theorem derived by Ger- 
schgorin [4] and rediscovered by Brauer [5] regarding 
the areas in the complex plane in which the roots of 
a matrix lie; it can be shown (see appendix, section 
3.3) that for our particular case, eq (5a), the eigen- 
values lie between zero and twice the most negative 
diagonal element in the matrix (a i3 ) derived from 
A in eq (5a) by setting D=0 and changing the sign. 
Since, by definition, kij>0, the roots are all negative. 

The proof that the eigenvalues are distinct if the 
kij are distinct has not been accomplished. 

From the theory of matrices we obtain two other 
properties of the characteristic determinant A. In 
the matrix {a i3 ), defined above, it is evident (see eq 
(5b)) that each column sum vanishes. This in- 
dicates that one eigenvalue of A corresponds to the 
root D—0. Next, the trace, or diagonal sum of a 
matrix (6 y ) is equal to the sum of the eigenvalues, 
XX> of the characteristic determinant \b i} — \SJ'|, 

s 

where 5} is the Kronecker delta. From eq (5) it is 
seen that for the matrix (—a tj ) the diagonal sum is 
yVflr™. where the prime means i^j. Since the eigen- 
O 
values of the characteristic determinant are un- 

8 .See appendix, section 3.2. 



changed if each element is replaced b} T its negative, 
the sum of the eigenvalues of A is — S'&«« The 

latter property is useful as a check when the eigen- 
values are determined in a particular problem. 

2.4. Relaxation Times 

With the general results obtained in section 2.3, 
we may write the solutions of the rate equations, 
noting that/i is zero, in the form 

Ni=Cn+C i2 eW+ . . . 

=Cn+C i2 e-^+ . . . (40 

where the rp are defined as the molecular relaxation 
times [1,2,3] 

re=-l/Mk). (6) 

Since we have shown that the fp are all real and 
negative, it follows that the rp for the general six-site 
Cartesian axes model are all real and positive. As 
implied earlier, eq (4') show how and at what rate 
a system tends to return to equilibrium after the 
abrupt removal of a disturbance such as a static 
electric field. Thus, for any arbitrary choice of the 
ley, consistent with eq (2), the return of the system 
to equilibrium is described by an exponential decay 
and cannot include any oscillatory modes. 

Also, from the application of Gerschgorin's Theo- 
rem to the range of the eigenvalues, we find that the 
shortest relaxation time r mln is 6 



Tmin = 



1 



: 2\a m 



(7) 



where a mm is the most negative diagonal element in 
the matrix (a i3 ) defined in section 2.3. If k pq is the 
largest of the k ih eq (5) and (7) lead to 



Train ^ 



N/'„ 



(8) 



The above results cover the inferences that may 
be readily drawn regarding the molecular relaxation 
times for the general six-site model. 

2.5. Dielectric Relaxation Spectrum 

If we assume the molecular and macroscopic 
relaxation times are equal, an approximation used 
previously [1,2,3,5], the molecular relaxation times rp 
describe the approach to equilibrium of the polariza- 
tion. The coefficients dp (for 0^2) in general 
measure the "strength" of each mode contributing 
to the relaxation process. In the case where we 
consider the restoration of equilibrium following the 
removal of a static electric field, the polarization 

6 In eq (7) the symbol |a mm | represents the absolute value of the a mm . 
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Pit) decays with time as 






- t/T e. 



(9) 



The coefficients Pp, and hence the orientational 
polarizabilities ap associated with each mode, may 
be obtained from the quantities (7^(/3^2) for a 
particular model by standard methods [1]. The ap 
are calculated assuming a polycrystalline substance. 
We note that for sinusoidal fields eq (9) leads to 
the approximate expression [1,6] 







(10) 



where e"(a>) is the dielectric-loss factor as a function 
of angular frequency w, K is a constant, and ap is 
the orientational polariz ability associated with the 
/3th mode. Although the expression is not exact, it is 
sufficient for an approximate analysis of the general 
characteristics of the dielectric relaxation spectrum. 
For the case of a single relaxation time, eq (10) 
reduces to the form of the Debye equation [7] for 

The coefficients dp (and hence the Pp and ap) 
cannot be obtained exnlicitly in the general case as 
they depend (eq (3) and (4)) on the eigenvalues, and 
the latter cannot be expressed explicitly in terms of 
the k tj . As a consequence, we have not been able to 
prove in the general case, where the k tj are distinct, 
that all relaxation modes are active in the dielectric 
spectrum. However, the latter situation seems 
quite possible, since in earlier calculations inactive 
modes (compensatory reorientation of dipoles to 
maintain a net dipole moment of zero) were formed 
only for very simple models, such as those containing 
two pairs of equivalent sites. 

To obtain further information about the dielectric 
behavior of dipolar systems represented by the gen- 
eral six-site model, it is necessary to examine special 
models. This is planned for a subsequent paper. 

A few useful formulas to apply to special models 
are given in the appendix, section 3.4. 

The author is indebted to John D. Hoffman for 
numerous helpful discussions regarding the general 
subject and to Karl Goldberg and Martin Greenspan 
for assistance with mathematical aspects. 

3. Appendix 

3.1. Relations Between k# 

Let the free energies relative to an arbitrary zero 
of energy for two adjacent sites i and j be designated 
as Vi and V h and the free-energy barriers for % to j 
and j to i transitions be defined as W tj and W jiy 
respectively. The height of the free-energy barrier 
between sites i and j, referred to the zero of energy, 
is then 

W{ij)=V t +W u =V,+W it . (11a) 

Consider a series of transitions i to j, j to I, I to i for 



a closed circuit of three adjacent sites and write the 
equations corresponding to eq (11a) for the barrier 
heights referred to the zero of energy 

j->l W(j,l) = V J +W jl =V l +W lj (lib) 

l->i W(Li) = V l +W li =V i +W il . (lie) 

The sum of the eq (11) is 

W v + W„+ W H = W jt + W tJ +W u . (12) 

If we recall the definition of k ih eq (1), inspection 
shows that eq (12) corresponds to log (£<;£; {&»)= 

log (kjikijku), and hence 



fc ij/Cj \k.n — tc i \k, ijfCj i . 



(13) 



The same type of relations, eq (12) and (13), can be 
obtained for any closed circuit of transitions between 
adjacent sites with three or more jumps to the cir- 
cuit. 7 For the six-site model there are seven inde- 
pendent relations, such as eq (13), namely, any 
seven of the eight eq (2). All others and those in- 
volving more sites are derivable from eq (13). Ac- 
cordingly, of the 24 transition probabilities only 17 
are independent. This is plausible on physical 
grounds. The local free-energy system of the general 
model is fixed when six minus one site energies and 
12 barriers W(i, j) are specified; the "minus one" 
arises because the zero of energy is arbitrary. 

The following geometrical argument illustrates 
how all product relations of the type in eq (13), but 
involving more sites, can be constructed from suit- 
able products of the product relations given in eq (2) 
for three adjacent sites. 

The six orientational sites (fig. 1) are located at the 
corners of an octahedron. Associated with the 
product relations on the k ih eq (13), for a closed cir- 
cuit of transitions between the adjacent sites i, j, Z, 
is the circuit of the perimeter of the face with vertices 
i, j, I, and the inverse circuit. Consider a closed 
circuit along the edges of the octahedron connecting 
the successively adjacent sites i, j, I, . . ., m, i, for 
example, where there are 4, 5, or 6 sites. This closed 
circuit divides the octahedron, a simply connected 
surface, into two parts, each of which is constituted 
of triangular faces; for convenience, choose the sur- 
face with the fewer faces, labeled T ly . . ., T p . If we 
traverse in succession in the same sense 8 as that on 
the whole perimeter (i, j, /,..., m, i), the perimeters 
of each of the triangular faces T 1} . . ., T v this is the 
same as traversing the bounding path i, j, I, . . ., 



7 An alternate method of arriving at product relations of the type in eq (13) or 
eq ^2) is to invoke "detailed balancing" for each pair of adjacent sites. This 
means that at equilibrium the number of dipoles turning in unit time from site 
i to site j is equal to the number turning from site./ to site i, or kaCn = kjiCji; On is 
the equilibrium number of dipoles in the ith site. For a particular case, sites 1, 2, 
3,wemaythenwritefci2C'ii = /c3iC2i,A;23C'2i=/:32C3i,and/:3iC'3i=A:i3C , ii. The product 
of the last three equations yields eq (2a) , fci2fc23fc3i=/c2i/c32/t - i3. 

8 That is, if the perimeter of the portion of the surface bounded by the sites 
(i, ;, . . ., m, i) is traversed in a clockwise sense by the above sequence, then each 
of the constituent triangles is traversed in the same sense. 
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m, i. Hence if we multiply the product equations 
such as (13) for the triangles T 1} . . ., T p , traversing 
their edges in the correct sense, we obtain the corre- 
sponding product relation for the circuit (i, j, I, . . ., 
m, i) and the inverse circuit (i, m, . . ., Z, j, i). In 
this way the k tj product relation for any sequence of 
connected sites can be formed by multiplying the 
product relations for the enclosed triangular circuits. 

3.2. Application of Goldberg's Theorem 9 to the 
Matrix for the Six-Site Model 

Theorem: If a real matrix A= (a tj ) of order n has 
the properties a ti >0 for i, j=l, . . . , n with i^j 
and 

fydftioii ' ' • a i m ii =a i4*Q / i-,io • • • ®i,i m (14) 



for all (i i} . . . , i m ) belonging to (1, . . . , 7i),m>3,then 
all the characteristic roots of A are real. 

We wish to apply the theorem to the matrix — A 
derived from A in eq (5a) by letting D be zero and 
changing the signs of the elements. If we let a tj be 
the element of — -Ao in the ith row and the^'th column, 
the usual order of index notation, we note that for 
is* j 

aij=kjij (15) 

i. e., the order of indices is reversed on the k's. Hence 
with eq (15), eq (14) can be written as follows, after 
interchanging the left and right sides and then re- 
versing the factor order on the right side: 

kiii 2 ki 2 h . . . ki m ii = ki 2 iiki- 3 i 2 • • . &Mm 

= ki 1 i m . . . ki z l 2 kl 2 i x . (140 

Equation (14') now evidently equates the product of 
ki/s for the closed circuit of m sites (#<ra<6) 
i\— H2, i2~ > is, • • • , im~^i\ to the corresponding 
product for the inverse circuit ii-*im, . . . , i2-n"i. 
There are two types of circuits: (a) circuits in which 
at least one successive pair i p , i Q are nonadjacent sites, 
and (b) circuits in which all successive site pairs are 
adjacent. Since by the S-J hypothesis kj \ = 
k iqip = for all nonadjacent site pairs, eq (14'), and 
hence eq (14), is valid for case (a). Next we have 
indicated in section 3.1 that closed-circuit product 
relations of the type eq (140 are valid for a sequence 
of three or more adjacent sites. Hence our matrix 
— Ao for the general six-site model satisfies the above 
theorem. 

3.3. Range of Eigenvalues 

Brauer [5] has derived the following theorem 
regarding the characteristic roots (eigenvalues) of a 
matrix. 



8 K. Goldberg of the Applied Mathematics Division of the Bureau. His proof 
is published in RP2652 in this Journal. The theorem as stated above is in a 
slightly different form from that as stated by Goldberg; his condition (2) has 
been replaced by the condition ai,->0 which still permits the construction of an 
equivalent symmetric matrix, namely, bn=(anaji)M. 



Theorem: Let A= (a f; ) be an arbitrary matrix and 

n n 

2-i\ a ij\—Pi) Z-j\ a ij\ = Qj' 
j=l i=l 

Each characteristic root X lies in at least one of the 
circles 

\z— a>u\<Pi 
and in at least one of the circles 

\z-djj\KQj. 

A is a square matrix of order n over the complex 
field; z is the complex variable x-\-iy. P t is the sum 
of the absolute values of the off-diagonal elements in 
the ith row; Qj is a similar sum for the jih column. 

In the matrix — A derived from the characteristic 
determinant, eq (5a), by letting D be zero and chang- 
ing the signs, we know that the off-diagonal elements 
are all real and positive, and so an off-diagonal 
column sum is Qj. Also from eq (5a) we see that 
the jth. diagonal element a^ is —Qj. Hence, for the 
jt\\ circle (j = l, . . . , 6) the radius is Q jy and the 
center is on the real axis at —Qj. Let Q m >Qj 
(j = l, . . . , 6). Q m defines the largest circle. Pre- 
viously all roots were shown to be real. Hence all 
roots lie on or within the mth circle on the axis of 
reals, that is between zero and 2a mm (a mm being the 
most negative of the diagonal elements in — A ). 

3.4. Useful Formulas 

Relation between the equilibrium number of dipoles 
in each site with and without a static electric field. We 
refer to the method of solving the rate equations in 
reference 1. After a "mode analysis" in which for 
each given mode the coefficients C& are determined 
relative to one another, the boundary condition at 
£=0 is introduced. This condition corresponds to 
the sudden removal of a static electric field. This 
results in the equations 

Cn = 2-jCi0, 

where Cfi is the equilibrium number of dipoles in the 
ith site in the presence of the field. A general 
expression for these coefficients in terms of the 
equilibrium coefficients for zero field is helpful in 
solving any particular model. 

The equations below for the Cf x are derived in a 
manner similar to that indicated in section 2.4 of 

reference 1. 

— * 

We assume a static electric field F defined by the 
magnitude F and the polar angles £ and 0, with £ 
being the azimuth angle measured from the 2-axis. 
For convenience, let 

y=/dF/kT and y x =y sin £ cos 6 
Jy—y sin £ sin 6 

y z = y COS f. 
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Then 



On — On 

ri f rt 

O21 — L/21 
641 =641 

ri f n 

^51 — ^51 
C/61 =C61 



l -7,(C 31 - C a )!N-y^C 2 i-C a )IN-y l ( 



^11 ^51 



JV 



0] 



1 -7* ( CSn ^ a - l)-T,(ai- <7«)/iV-%(C„-C 61 )/iv] 
l-7^C„-Ck 1 )/A r -T,(^^i+l)-Y.(Cit-Q B )/iv] 

1-7, (^ i ^+l)-7„(ai-<7 tt )/iV-7 2 (a i -<7 M )/iv] 

1 -7,«7 31 - Cn)/N~y v (C 21 -C tt )IN-y t ( gll ~ (761 +l)] 



TAe fofoZ polarizability, a tot . We assume a poly- 
crystalline substance and follow the procedure out- 
lined in section 2.4, reference 1, for obtaining the 
polarizability associated with each mode. Here, 
however, we introduce for each site i the value of 
Ni at t=0, namely, CJi. With the use of the ex- 
pressions for the Cfi given in the preceding para- 
graph, one can obtain 



Oitot=\ 1 



(g31-g»l) 2 +«721"g4l) 2 +(Cll~g«lH M 1 



iV 2 



] 



3kT 

(17) 



(16) 
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